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Using Monte Carlo integration techniques, we investigate running coupling 
effects compatible with the high energy bootstrap condition to all orders in 
the strong coupling in evolution equations valid at small values of Bjorken x 



in deep inelastic scattering. A model for the running of the coupling with an- 



alytic behavior in the infrared region and compatible with power corrections 



to jet observables is used. As a difference to the fixed coupling case, where the 



momentum transfer acts as an effective strong cut-off of the diffusion to in- 
frared scales, in our running coupling study the dependence on the momentum 
transfer is much milder. 



1 Theoretical set up 

In this work we present a study of a long-standing problem in high en- 
ergy QCD: the treatment of the running of the strong coupling in evolu- 
tion equations driven by the Balitsky-Fadin-Kuraev-Lipatov (BFKL) equa- 
tion [TH3]. This subject has been extensively discussed in the literature (see, 
e.g. Refs. [1]). We put the emphasis on a particular form of writing the equa- 
tion with running coupling which is consistent, in principle to all orders in a 
coupling expansion, with the bootstrap property of QCD scattering amplitudes 
at high energies [5HZ] - Bootstrap in high energy QCD has been discussed in 
many papers [8] and we refer the reader to the original work by Braun [5] for a 
detailed discussion directly related to our present study. Our new contribution 
is to be able to study the problem using Monte Carlo integration techniques 
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which solve the BFKL equation with a running coupling exactly and allow us 
to access exclusive information of the final states in the cut amplitude case, 
and of the diffusion pattern in the virtual diagrams for the non-forward elastic 
amplitude. We are particularly interested in the dependence of our solution 
on the total momentum transfer. For a connection of this representation of 
the BFKL equation with renormalon contributions we refer the reader to [5] 
and for a more recent related analysis in coordinate space to [7]. 

To model the running of the coupling in the infrared we will use a simple 
parametrization introduced by Webber in Ref. [9] which at low momentum 
scales is consistent with global data of infrared power corrections to pertur- 
bative observables (mainly related to jet event shapes). The relevant formula 
reads 



47T 

a s (k) = — , 
1 ' A> ling 



1 125 (A 2 + Ak 2 



(A* - P) (4 + g) 



which, for (3 = (llN c -2n f ) /3, n f = 3 and A = 0.25 GeV, gives a s (91GeV) 
0.118. Its dependence on k is shown in Fig. [TJ 



Model for the running coupling in the infrared 




k (GeV) 

Figure 1: Model for the running of the coupling in the infrared. 



In the following we introduce the parametrization ([I]) of the running of the 
coupling in a version of the BFKL equation for a projection on a SU(3) singlet 



2 



in the t-channel which can be directly derived from Eq. (4.12) in Ref. [luj^l . If 

d 2 k Q 2 



we US' 



a. 



A; 2 (A; — q)' 



(2) 



as the notation for the gluon Regge trajectory, we can then write the non- 
forward BFKL equation with a fixed coupling for the t-channel partial wave 
fu(fa, fa] q) as 



cj-e{fa)-e{fa-q) )Ufa,fa;q) = 5^ 2 \fa - fa) 



(I ~ q) 2 kf (I - faf 



(fa - qfl 2 
(I - qfk 2 



1 + 



f u (i,fa;q), ( 3 ) 



with a s = a s N c /n. Note that the solution to this equation corresponds to 
a four point Green function for four off-shell reggeized gluons carrying two- 
dimensional transverse momenta —fa, fa — q, fa, q — fa, all of them outgoing. 
q corresponds to the total momentum transfer in the t-channel. In order to 
match the normalization used by Braun in [5] it is now needed to introduce 
the rescaling 



I 2 

G u (l,k 2 ;q) = fu(l,k2]q)jp 



(4) 



The set up under study can be simply implemented by replacing each of the 
squared transverse momenta p 2 of Eq. (E]) with a general function rj(p). The 
new trajectory will then be 



<q) 



d 2 k 



77(g) 



Atx r](k)rj(k — q) 



(5) 



and Eq. (j3J) now reads 

to - e{fa) - e{fa - q) ) G w {fa, fa; q) = 5^ (fa - fa) 



+ 



d 2 l r)(fa) 
2n T](l)r](fa - I) 



v(fa - gMQ - v(q)v(fa - 1) 

rj(l-q)rj(fa) 



g u {l,fa;q)- (6) 



1 For a Monte Carlo study of the total momentum transfer dependence of the BFKL 
gluon Green function in the LO and NLO adjoint representations in QCD and N = 4 SUSY 
see Refs. [11] 

2 In this text all momenta are two-dimensional and we take the notation p = p, for moduli 
we use \p\ 
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As it has been shown in [5J, this equation is compatible with the bootstrap 
condition for the all-orders expansion of the function rj. 

Let us first study the simplest case, with q = 0, corresponding to forward 
scattering, or, by the optical theorem, to a contribution to the total cross 
section. From now on we also fix r)(k) = k 2 /a s {k) with a s as in Eq. (pQ). We, 
therefore, can write 

u -2e(fa)^G UJ (k 1 ,k 2 ) = 5 (2 \fa - k 2 ) + J ^ g u {l,k 2 ) (7) 

since 77(0) = 0. 

For our Monte Carlo implementation of the solution to this equation (see 
Refs. [12] for similar studies in the fixed coupling case) it is convenient to 
introduce a shift in the integration momentum of the form I = k + fa and a 
mass parameter A to separate the resolved real emissions (with k 2 > A 2 ) from 
the unresolved ones (with k 2 < A 2 ). The latter, after integration over the 
phase space of the emitted gluons, generate infrared divergences which should 
cancel against those of the gluon Regge trajectories. The final results here 
presented are independent of A in the limit A — > 0. Taking the approximation 
Gu{k + ki, k 2 ) — Gw{fa, fa) for unresolved emissions we can then write 



(^-2e x {fa)^GUfa,fa) = ^\fa-k 2 ) 



tPk r,(fa)9(k?-\*) n 



where 



_ fd 2 k V (q)6(k 2 -\ 2 ) 
X[q) J 2nr ] (k)(r ] (k)+ V (k-q)Y 1 1 

We go back to Bjorken x space using 

F(fa,k 2 ,x) = J ^L x -"Gu{ki,k 2 ). (10) 

The final expression to be evaluated using Monte Carlo integration techniques 
is 

Pi 



{ ~ jl r d 2 r>- v (fa + Ez-i 1 

r(fa, k 2 , x) = x- 2 ^ (fa -k 2 )+j2U — — v 7 : ^ 

l n=Xi=l J * v(Pi)V (fa + Yd=lPl^ 

x6(p 2 -\ 2 )5v(fa + ± Pl -k 2 ) f ^^(^^^^(^SU^JLuj 
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where Xq = x. Note that n corresponds to the number of on-shell gluons 
emitted with a longitudinal momentum fraction Xi and a transverse momentum 
k% ■ 

For completeness, we compare the gluon trajectory in the form 



e A(<?) 



q 2 f d 2 k a s (k)a s (k — q) 



a s (q) J 2itk 2 k 2 + (k-q) 2 



9(k 2 -X 2 ), 



(12) 



with the usual one at leading order in our A-regularization of infrared diver- 
gencies 



2nk 2 k 2 + (k- q) 2 



«»(g) n 9 2 hQ x 
^~ ln p (13) 



in Fig. |2] (both lines are calculated with A = 0.01 GeV). Note that the be- 
haviour of both representations is quite different at large values of the mod- 
ulus of the transverse momentum in the reggeized gluon propagators (which 
correspond to the power-like terms in Eq. (111))). 




q (GeV) 

Figure 2: Gluon Regge trajectory in different schemes 
A similar iteration of the BFKL kernel applies for the more complicated 
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non-forward equation. In this case Eq. ( II ip should be 



J r (k 1 ,k 2 ,q,x) 



X 
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(14) 
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X i 
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The regularization of infrared divergencies is identical to the forward scattering 
case, with a unique infrared regulator A. Clearly the \q\ — > limit of Eq. (fT5|) 
corresponds to Eq. (TTTT) . We are now ready to present the numerical results 
for the evaluation of both functions in Eqs. (TTTj) and (|T5|) . 

2 Numerical results & discussion 

We start by investigating the evolution with x of the gluon Green function 
for fixed coupling and show our results in Fig. [3j In the region of x where we 
have performed our calculations (x > 0.0003) we observe that the growth of 
the function as x decreases is more pronounced for small q. Due to SL(2, C) 
invariance the asymptotic slope of our curves will be the same in the x — Y 
limit. We have presented the plot for a s = 0.2, \ki\ = 10 GeV and |&2| = 
15 GeV, but our results are generic. 

As it is well-known, the effect of introducing the running of the coupling is 
to reduce the growth of the Green function as x goes to zero. The advantage 
of our method of calculation to previous analysis in the literature is that we 
can solve the BFKL equation exactly, with no asymptotic approximations. 
In general, we find that the introduction of the running of the coupling in a 
form compatible with bootstrap to all orders indeed reduces the growth of the 
solution. As in the fixed coupling case, with running we see in Fig. 2] that we 
have a smaller Green function as the momentum transfer q grows. However, 
in the running coupling case we have to go to smaller values of x if we want 
to make this feature manifest. 
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Figure 4: Evolution of the gluon Green function with the Bjorken variable x, 
for fixed values of the transverse momenta and a running coupling. 
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In order to investigate this effect in more detail we can perform an analysis 
of the gluon Green function in the collinear regions where \ki\ 3> \k 2 \ and vice 
verse. Let us fix \k 2 \ = 20 GeV, and study the solution to our equation for 
fixed and running couplings, with different values of x. In Fig. [5j for a fixed 
coupling and a large value of x (x = 0.135) it can be seen that the effect of 
increasing the value of q is important only in the region \ki\ <C \k 2 \. However, 
when we go to higher energies and decrease x, as in Fig. El we see that the 
effect of q is manifest in the full range of \k\\. 

Fixed coupling a = 0.2, q = GeV 

q = 3 GeV 

q = 6 GeV 



10 15 20 25 30 35 40 

Figure 5: Collinear behaviour of the gluon Green function, for fixed values of 
one transverse momentum, a non-running coupling and x = 0.135. 

It is striking that the running of the coupling completely eliminates any 
effect of introducing the momentum transfer for x = 0.135 and q from zero 
up to 6 GeV. This can be seen in Fig. [3 It is needed to reduce the value of 
x by one order of magnitude in order to start finding a decrease of the Green 
function as q increases, for < \k 2 \, as in the fixed coupling case. This is 
shown in Fig. [HJ 

Since we have full access to the exclusive information of all momenta 
configurations in the gluon ladder let us investigate the typical transverse 
momentum scale running in the internal propagators. This corresponds to 
the well-known diffusion "cigar" [13] where the mean value of the variable 
r = log < pf > / (GeV 2 ) is plotted (together with the lines of one standard de- 
viation towards the infrared and ultraviolet for the gluons in the BFKL ladder) 
as a function of the normalized rapidities of the corresponding gluon lines. For 
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Figure 6: Collinear behaviour of the gluon Green function, for fixed values of 
one transverse momentum, a non-running coupling and x = 0.018. 
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Figure 7: Collinear behaviour of the gluon Green function, for fixed values of 
one transverse momentum, a running coupling and x = 0.135. 



9 




Figure 8: Collinear behaviour of the gluon Green function, for fixed values of 
one transverse momentum, a running coupling and x = 0.018. 



a s = 0.2 in Fig. |9]we show the typical effect of introducing a non-zero momen- 
tum transfer: the diffusion to the infrared is reduced as q increases, while the 
ultraviolet diffusion line remains stable. In other words, the momentum trans- 
fer acts as an effective infrared cut-off. For the same large value of x — 0.37 
we find in Fig. [10] a smaller suppression of the infrared diffusion in the set up 
with a running coupling as the momentum transfer increases. We also see that 
the diffusion to the ultraviolet is suppressed with respect to the fixed coupling 
case, independently of the value of q. 

For a smaller value of x the spread in transverse momentum in the internal 
gluon propagators is much bigger, in particular for a fixed coupling, see Fig. [TT1 
The modification of the diffusion picture depending on the values of q is inter- 
esting since the suppression of the evolution into the infrared is very big but 
it is also present with the opposite effect, an enhancement, in the ultraviolet 
region. In the latter a large momentum transfer pushes the gluon momenta to 
live in more perturbative regions of phase space. This diffusion to the ultra- 
violet is somewhat reduced in the case with a running coupling, see Fig. [12j 
As the value of x gets smaller we can see that the influence of introducing a 
non-zero momentum transfer is larger, always "pulling" the diffusion "cigar" 
towards more perturbative regions in both, the fixed and running coupling 
scenarios. 
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Figure 9: Distribution of the transverse momenta in the internal propagator 
of the gluon ladder for a fixed coupling and a large value of Bjorken x. 




Figure 10: Distribution of the transverse momenta in the internal propagator 
of the gluon ladder for a running coupling and a large value of Bjorken x. 
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log (1/x)/log (1/0.05) 

Figure 11: Distribution of the transverse momenta in the internal propagator 
of the gluon ladder for a fixed coupling and a small value of Bjorken x. 




Figure 12: Distribution of the transverse momenta in the internal propagator 
of the gluon ladder for a running coupling and a small value of Bjorken x. 
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3 Conclusions & outlook 

We have investigated small x evolution equations which incorporate a running 
coupling compatible with bootstrap to all orders in a perturbative expansion. 
Our analysis is exact since we have found the solution to the BFKL equation 
using an iteration in transverse momentum space which is finally expressed 
in terms of integrals over transverse momenta and the rapidity of the internal 
gluon propagators, which we have evaluated using Monte Carlo integration 
techniques. We have found that the effect of the total momentum transfer 
as an effective cut-off for diffusion of transverse momenta into the infrared 
region is very suppressed in this way of running the strong coupling. The 
propagation of the gluon ladders into ultraviolet regions is also suppressed, 
independently of the value of the momentum transfer. Our next task will be 
to integrate these results with suitable impact factors in order to gauge the 
phenomenological relevance of the results here presented. 
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